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14 Abstract 

15 We report analyses of a dataset spanning 38 years of near-annual fishing for Antarctic toothfish 

16 Dissostichus mawsoni, using a vertical setline through the fast ice of McMurdo Sound, 

17 Antarctica, 1972-2010. This constitutes one of the longest biological time series in the Southern 

18 Ocean, and certainly the longest for any fish. Fish total length, condition and catch per unit 

19 effort (CPUE) were derived from the >5500 fish caught. Contrary to expectation, length- 

20 frequency was dominated by fish in the upper half of the industrial catch. The discrepancy may 

21 be due to biases in the sampling capabilities of vertical (this study) versus benthic (horizontal) 
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22 fishing gear (industry long lines), related to the fact that only large Antarctic toothfish (>100 cm 

23 TL) are neutrally buoyant and occur in the water column. Fish length and condition increased 

24 from the early 1970's to the early 1990s and then decreased, related to sea ice cover, with lags 

25 of 8 months to 5 years, and may ultimately be related to the fishery (which targets large fish) 

26 and changes in the Southern Annular Mode through effects on toothfish' main prey, Antarctic 

27 silverfish Pleuragramma antarcticum. CPUE was constant through 2001 and then decreased 

28 dramatically, likely related to the industrial fishery, which began in 1996 and which 

29 concentrates effort over the Ross Sea slope, where tagged McMurdo fish have been found. Due 

30 to limited prey choices and, therefore, close coupling among mesopredators of the Ross Sea, 

31 Antarctic toothfish included, the fishery may be altering the trophic structure of the Ross Sea. 
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51 

52 Introduction 

53 Unlike any other marine ecosystem in the world, the high latitude waters of the Southern 

54 Ocean contain a single lineage of teleost fishes that dominates fish diversity, abundance and 

55 biomass (Eastman 2005). Tectonic, oceanographic and climactic changes over the past few tens 

56 of millions of years eliminated previous fish faunas and provided an opportunity for this lineage 

57 of perciform notothenioids to radiate opportunistically and to fill most benthic and pelagic 

58 niches, despite the lack of a swim bladder. For example, the 2- meter-long Antarctic toothfish 

59 Dissostichus mawsoni, and the 15-25 cm-long forage species, Antarctic silverfish Pleuragramma 

60 antarcticum, are key piscine predator and prey in the ecosystem in spite of the incongruity of 

61 being closely related sister species/genera (Balushkin 2000). As adults Antarctic toothfish are an 

62 important component of the community of mesopredators (also including penguins, seals, and 

63 cetaceans) that inhabit the Ross Sea, Southern Ocean; and most of them, too, depend on 

64 silverfish (Ballard et al. 2011). To date, the Ross Sea is the least anthropogenically altered 

65 stretch of ocean remaining on Earth (Halpern et al. 2008) and, thus, the interactions of this 

66 group of mesopredators and their prey are of particular interest as a model for how cold-water, 

67 continental shelf food webs once operated elsewhere. 

68 In this system, productivity is highest and trophic interactions most intense during 

69 summer. Most of the mesopredators are present over the Ross Sea (shelf and slope) only 

70 during summer, moving north with daylight and the pack ice during winter; a few —Antarctic 

71 toothfish, emperor penguin Aptenodytes forsteri and Weddell seal Leptonychotes weddellii — 

72 remain year round. All feed during summer primarily on three prey species over the shelf: 
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73 Antarctic silverfish and two species of krill: Antarctic krill Euphausia superbo and crystal krill E. 

74 crystallorophias (Ballard et al. 2011, Smith et al. 2011). Over the southern shelf, the silverfish by 

75 its numbers and diet is thought to be the main predator of crystal krill and the toothfish to be 

76 the main predator of silverfish— as measured in McMurdo Sound, 71.2% by frequency of 

77 occurrence and 89.2 % by dry weight (Eastman 1985b, La Mesa et al. 2004), consistent with the 

78 generalization that fish are the most important predators in most marine ecosystems (Sheffer 

79 et al. 2005). Off the shelf, the diet of these mesopredators is much more diverse (Ainley et al. 

80 1984, Fenaughty et al. 2003). Due to the reduced number of dietary items over the shelf, a tight 

81 trophic coupling has been observed among them, e.g. prey consumption of cetaceans 

82 negatively affects the foraging area and diet of penguins (Ainley et al. 2006), and the foraging of 

83 Weddell seals can negatively affect the prevalence of toothfish (Testa et al. 1985), although the 

84 consequences of this fish-seal interaction for local availability of silverfish remain unknown. 

85 In this system, too, the apex predators are the Ross Sea killer whale Orcinus [ orca ]; 

86 thought to be a separate but yet-to-be-described species; Morin et al. 2010), which has been 

87 observed to feed substantially on toothfish (Ainley et al. 2009); and the larger killer whale 

88 "ecotype B" (O. orca), which feeds on seals and likely large penguins (Pitman and Ensor 2003, 

89 Ainley and Ballard 2011). Weddell seals also prey to a significant degree on toothfish (Ainley 

90 and Siniff 2009, and references therein; Kim et al. 2011), and along the Ross Sea shelf break and 

91 north, toothfish are fed upon by southern elephant seals Mirounga leonina, sperm whales 

92 Physeter macrocephalus and colossal squid Mesonychoteuthis hamiltoni (Yukhov 1970, 

93 Pinkerton et al. 2010, Ainley 2010, Smith et al. 2011). 
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94 To understand the trophodynamics of this system on the Ross Sea shelf requires 

95 detailed information on all of its predators and mesopredators. While much is known about the 

96 air-breathing members of the group, all of which have nowhere to hide from researchers, and 

97 indeed have been well researched in the Ross Sea since the mid-1900's, learning about 

98 toothfish life history and abundance has been a challenge. On the one hand, the morphology 

99 and physiology of Antarctic toothfish are among the best known of Southern Ocean species: it 

100 is one of just five neutrally buoyant nototheniids (Eastman and DeVries 1981), attaining this 

101 status at around 100 cm TL (Near et al. 2003), and ecologically dominates the Southern Ocean 

102 fish fauna in the sense that it is the major piscine predator (Eastman 1993). It has a suite of 

103 adaptations to allow presence at subzero temperatures including blood antifreeze (DeVries 

104 1988, DeVries and Cheng 2005) and other adaptations (Eastman 1993); it grows rapidly when a 

105 subadult, then its growth rate slows, and it can live to 50 yrs. On the other hand, the ecology 

106 and population dynamics of Antarctic toothfish have remained obscure. From a fishery, the 

107 "legal" portion of which began in 1996-97 (hereafter we will identify austral summers by the 

108 initial year, 1996 in this case), we have learned about the geographic aspects of size-frequency 

109 distribution (Hanchet et al. 2008, 2010), confirmed growth rates (cf. DeVries and Eastman 1998, 

110 Horn 2002, Brooks et al. 2010) and gained insights into variation in condition (Fenaughty et al. 

111 2008), age of recruitment (Parker and Grimes 2010) and diet of individuals in the deepest parts 

112 of the Ross Sea along its continental slope (Fenaughty et al. 2003). Otherwise, much remains to 

113 be learned about the Antarctic toothfish, unlike its only congener, the non-Antarctic Patagonian 

114 toothfish (D. eleginoides), about which much has been learned (Collins et al. 2010). Antarctic 

115 toothfish are thought to spawn during winter; no free eggs or larvae and rarely fish <50 cm 
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116 have been collected in the Ross Sea region and the spawning frequency, fecundity and aspect 

117 of larval and juvenile life history are unknown. This information gap is the result of severe 

118 restrictions on scientific investigation resulting from this Antarctic species' "preferred habitat": 

119 deep, cold, sea ice covered ocean. Any data from the fishery are limited to just the 3-4 months 

120 of ice-free summer (Dec-March). 

121 Toward learning more about this important Southern Ocean mesopredator, we present 

122 results of analysis on data collected from >5500 fish captured, measured, marked and released, 

123 during an extended effort led by A. DeVries almost annually in McMurdo Sound, southern Ross 

124 Sea, 1972-2010. The data set thus constitutes one of the longest biological time series available 

125 for the Southern Ocean, the only long one for a fish species, and one that persisted through a 

126 period that has seen interesting changes in the Ross Sea, in large part driven by climatic forces. 

127 The Southern Annular Mode midway in the period (i.e, about 1986) shifted to mostly positive 

128 leading to a marked increase in winds, sea ice extent and persistence of the coastal polynyas 

129 that strongly affect Ross Sea processes (Parkinson 2002, Zwally et al. 2002, Jacobs et al. 2002, 

130 Stammerjohn et al. 2008, Jacobs 2006, Russell et al. 2006). Several of the air breathing 

131 mesopredators are known to have responded through changes mediated by sea ice (Ainley et 

132 al. 2005). In addition, minke whales Balaenoptera bonaerensis were heavily hunted in the late 

133 1970s-early 1980s, affecting penguin competitors, but since have recovered; Weddell seals in 

134 McMurdo Sound, too, were also severely hunted during that period but also appear to have 

135 mostly recovered (Ainley et al. 2007, Ainley 2010). In 1996, a toothfish fishery was inaugurated 

136 in the Ross Sea region, including to within ~60 km of McMurdo Sound (opposite side of Ross 

137 Island; FAO/CCAMLR Area 88), and this is now the largest fishery for toothfish south of the 
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138 Antarctic Polar Front (CCAMLR 2010). Therefore, in our analysis, given the close coupling 

139 observed among air breathing mesopredators, we hypothesized that we should see trends in 

140 toothfish condition, size and abundance that correlate with the above described environmental 

141 and anthropogenic changes. After examining trends over time, we explore the significance of 

142 environmental variables, especially those related to ice cover, to gain insight into the extent to 

143 which these factors could be explaining the observed temporal patterns in toothfish 

144 characteristics. 

145 Methods 

146 Fish capture 

147 Scientific fishing occurred primarily from one site, about 4 km west of McMurdo Station, in the 

148 vicinity of 77° 51' S, 166 0 40' E (Fig. 1); it was accessible from the Station by a ~10 min over-ice 

149 drive by tracked or wheeled vehicle. Note that most of the study period occurred before the 

150 age of GPS and the ability to precisely identify the location of the fish hut was limited. We know 

151 that the fishing site shifted slightly from year to year, e.g. to deal with pressure ridges, snow 

152 banks etc, and after 1999, to accommodate aircraft approach to the Sea Ice Runway. Both the 

153 runway and then the fish hut were repositioned anew each season, as noted. In 2000, the 

154 runway was shifted more than usual, this time northward, to be closer to and in view of 

155 McMurdo Station, forcing movement of the fishing site more west and south than usual. 

156 Nevertheless, all sites but one were within about a 2 km radius, with depth ranging 415 - 495 

157 m. Moreover, all sites were within the area in which the survey by Testa et ol. (1985) logged 

158 high catch rates of Antarctic toothfish and the area within which Weddell seals have been 
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159 frequently seen with toothfish (Ainley and Siniff 2009), even to 2010. Beginning in 2008, the 

160 fishing was based out of Scott Base instead of McMurdo, requiring a longer drive by 

161 researchers. 

162 Using a small hydro-winch with 3/32-inch wire and a 25 kg weight attached to the 

163 bottom end to keep the wire under tension (i.e., a vertical set line), a line was deployed through 

164 a hole drilled through the annual fast ice. A heated fish house was placed over the hole. Pulling 

165 the line, depending on number of fish caught, took 2-3 h. Initially 15-22 one-meter leaders with 

166 swiveled, stainless steel, long shank #10 hooks spaced 20-25 m apart were fished (thus 

167 sampling the lower ~300 m of water column, but not waters within 10-20 m of the bottom to 

168 avoid scavenging by benthic amphipods). Quickly it was found that most of the fish were caught 

169 in the lower half of the hook array. Thus leaders were shortened to 30 cm and hooks were 

170 spaced 3-5 m apart, starting 10 m from the bottom (thus the lower ~100m of water column was 

171 sampled but, again, not the bottom). For the first 3 years, live Pagothenia borchgrevinki were 

172 used for bait, but it was then found that dead bait worked just as well and so thereafter New 

173 Zealand yellow eyed mull et Aldrichetta forsteri, cut in half, was used. Sometimes the interval 

174 between deployments ran 12 h but the interval for the very large majority was 24 h sets. Some 

175 ran 48 h if poor weather prevented access to the fish house, but the fish caught were often 

176 exhausted after such prolonged time on hooks. Once this was realized, this long a soak time 

177 was avoided if possible. When the fishing became based out of Scott Base, this necessitated a 

178 much longer commute to the fish hut, and these sets often soaked for >45 h. 
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179 Captured fish not used for experiments were placed in a V-trough (with a seawater 

180 soaked cloth over their eyes), measured to the nearest cm; weighed to the nearest pound 

181 (converted to kg); tagged with a numbered "Floy" dart tag behind the 2 nd ray of the 2nd dorsal 

182 fin and a tail locking tag; in many cases injected with tetracycline; and released. Tetracycline 

183 served as an otolith annuli marker in case the fish was recaptured. This procedure, lasting 3-5 

184 min, was done on the floor of the heated fish hut so that the fish neither warmed nor froze. The 

185 open sea surface was 0.5 m below the level of the floor so that the fish could be gently lifted by 

186 their gill covers and returned to the water without abrading their skin and causing scale loss. 

187 Small fish on the line for any length of time appeared to have been prone to attack by 

188 large fish, as those retrieved often had many parallel teeth marks on both sides of their body as 

189 well as many missing scales. This is consistent with reports of small toothfish found in the 

190 stomachs of larger ones taken in the industrial fishery (Petrov and Tartarnikov 2010). Because 

191 standard, stainless steel, long-shank hooks were used (not short-shank with the gap of the tip 

192 being less than that at the lower part of the curve, as in the industrial fishery; CCAMLR 2008), 

193 most of the fish >160 cm were likely lost, given that some hooks were straightened out or 

194 broken. Although Weddell seals cannot entirely be ruled out as taking the bait or hooked fish, it 

195 certainly cannot be other fish species because there are none that are anywhere near the size 

196 of the toothfish in McMurdo Sound. When we retrieved toothfish in the presence of Weddell 

197 seals sharing the fishing hole, the seals paid little attention, apparently not recognizing the 

198 toothfish as prey in that context. [NOTE: We were always sufficiently far from the ice edge that 

199 killer whales would not have access to the fishing holes and we did not fish when the killer 

200 whales were present in the ship's channel.] On the few occasions when a large treble or 
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201 industrial-type hook was used, more of the larger specimens were caught. These hooks, 

202 however, left large wounds in the jaw likely leading to reduced survival. Therefore, their use 

203 was discontinued. 

204 

205 Data analysis 

206 Fish metrics 

207 Overall strategy. Our primary objectives with regard to analysis were two-fold: first to 

208 characterize change in fish metrics (length, condition, and capture rates) over the study period 

209 (1972-2010), and then to identify environmental variables influencing variation in these metrics 

210 among years and across decades. Analyses were conducted on trends in condition and length 

211 of toothfish (for condition, n = 5,403 individuals, and for length, n = 5,437, among the 5,587 

212 caught among 27 yr of data), using both (1) mean condition index per year and (2) 50 th , 75 th , 

213 and 95 th percentiles of total length (TL) in each year. (3) We confirmed the analysis of annual 

214 values of condition and length by analyzing trends with respect to the condition and length of 

215 individual fish. Finally (4), we assessed number of individuals caught per month and year as part 

216 of a "catch per unit effort" study (see "Catch per unit effort" below). Fish having partial data or 

217 measurements that clearly were wrongly recorded were excluded from the data. 

218 We present results of toothfish length and condition in relation to variation among 

219 years, examining both linear and non-linear trends. To allow for non-linearity, we fit two types 

220 of models: (1) polynomials of second-order, and (2) results of "change-point" analyses using 


221 


linear splines (Harrell 2001, pp. 18-19). In the latter case, two linear segments (each with its 
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222 own slope) are joined at a specified "knot". This allowed us to estimate and test for changes in 

223 (linear) trend. Both the second-order polynomial and the "two-joined-linear-segments" model 

224 have two degrees of freedom and thus are easily compared with respect to fit (e.g., by 

225 minimization of deviance). We also explored other simple models, such as cubic polynomial and 

226 "step-function" models (where there is a disjunct step-wise increase or decrease). The favored 

227 function fit to the data, and presented here, was the one that optimized AIC (Akaike 

228 Information Criterion). 

229 Because the toothfish appears to prey primarily on silverfish over the shelf (see 

230 Introduction), especially those of approximately 4-5 years of age (Eastman, pers. obs., see 

231 below), and because La Mesa et at. (2010) noted an effect of the amount of open water on the 

232 production and survival of silverfish larvae, we considered environmental variables with lagged 

233 effects of 1 year (specifically, 8 mo) to 5 years. Therefore, we looked at the effect of ice cover 

234 and its persistence previous to the current (ice-covered) fishing season. 

235 

236 Condition. We calculated a Fulton-type condition factor (Anderson and Gutreuter 1983), 

237 scaled to center around 1.0, as an index of weight per unit length, an indirect estimate of fish 

238 girth (Davidson and Marshall 2010) and a proxy for body shape: 

239 K=(W/TL s )x 10 2 

240 where W = body weight in g and TL = total length in cm. Larger values for K indicate greater 

241 weight per unit length associated with a thicker body. We did not employ K as an indicator of 

242 body fat content. The resulting index varied from 0.577 to 2.99, with mean = 1.265 ± 0.149 SD. 
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243 Fenaughty et al. (2008) identified an important condition threshold, K ah , which designates "axe- 

244 handle" fish that are in particularly poor physiological condition, K ah = 1.01248. Therefore, in 

245 addition to changes in mean condition, we analyzed the proportion of individuals that each year 

246 was below this threshold. The term is derived from the long-thin body with a large head of 

247 emaciated fish. Higher K is generally associated with better physiological condition in wild fish. 

248 Individual condition data came almost entirely (99.9%) from captures in four months: 

249 September, October, November, and December; 28.1% of captures were in October, 50.0% 

250 were in November. Since captures were essentially confined to four months, we have analyzed 

251 monthly variation in condition and length by treating month as a factor with four levels (i.e., as 

252 a categorical variable). Thus, in controlling for the effect of month, we make no assumptions 

253 about how the dependent variables change from month to month. The sample size of captures 

254 per year differed greatly (Table 1). 

255 Mean condition index, K, was calculated for each year, as was the proportion offish 

256 below the axe-handle threshold. Total length was considered to be both an index of body size 

257 and a proxy of approximate age (Horn 2002, Brooks et al. 2010). To assess changes in the 

258 frequency distribution of fish length, with specific attention to fish that were of greatest length 

259 (and thus of greater age), we analyzed three metrics: the 50 th percentile (i.e., median), 75 th 

260 percentile (i.e., upper quartile), and 95 th percentile. The latter was chosen as a more statistically 

261 robust measure of "especially large" individuals than the simple maximum; the 95 th percentile 

262 reflects a larger sample (in fact, in 20 of 27 yr, at least 5 individuals were >95 th percentile). 
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263 Variation in the year-specific percentile values was compared to the percentile values for the 

264 entire study. 

265 For analysis of annual change in condition, we weighted the mean value by the inverse 

266 of the standard error of mean condition. Similarly, for annual change in length, we weighted 

267 the 50 th , 75 th , and 95 th percentiles for a given year by the inverse of the standard error of mean 

268 length for that year. Thus the weighting reflected the sample size as well as the individual-level 

269 variability in the metric for that year (e.g., the less variable, the greater the weight). 

270 Catch per unit effort. To analyze patterns in toothfish catch rate, i.e., catch per unit 

271 effort (CPUE), we analyzed the number offish caught divided by the most appropriate measure 

272 of capture effort. Among effort metrics evaluated, the metric most closely reflecting the 

273 number offish caught per year was the total number of hooks set in each year. Comparing 

274 results among years, the number of hook-hours did not correlate with the numbers of fish 

275 caught, once we controlled for the number of hooks set per year. We also analyzed effort as a 

276 separate covariate, where the dependent variable was the number of individuals caught. 

277 Results were similar using either method (dividing number caught by effort or statistically 

278 controlling for effort); we only report the former because it is more widely used. 

279 We conducted the abundance analysis in two ways: (1) captures per 10 hooks, natural 

280 log-transformed, and (2) proportional change in abundance from one year to the next. Note 

281 that the proportional change in abundance from year t, N t , to abundance in year t+1, N t +i, is 

282 equal to the antilog of [ln(N t +i )-l n(N t )] . Hence, we present graphically the results of capture per 

283 10 hooks in terms of In-transformed capture rates. 
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284 Effort differed among months, and so did the capture rates with respect to month (see 

285 below; also Testa et al. 1985). We therefore analyzed between-year variation in captures per 

286 unit effort while controlling for month of capture, and present results that reflect the statistical 

287 adjustment with respect to month. 

288 Furthermore, month of capture could confound observed patterns with respect to 

289 among-year variation in condition and length, given that monthly effort varied among years. To 

290 guard against such possible confounding, we repeated the analyses of condition and length, 

291 using individual observations instead, in relation to between-year patterns and with respect to 

292 the environmental variables. In all cases the pattern observed using annual measures of 

293 condition and length was confirmed as significant by the analyses of individuals, after 

294 controlling for month of capture. 

295 

296 Environmental variables 

297 Sea Ice. Using satellite passive-microwave data for the Ross Sea sector from NASA 

298 Nimbus 7 satellite and the Defense Meteorological Satellite Program (e.g., Parkinson 2002, 

299 Zwally et al. 2002), ice area and ice extent were calculated daily or every other day for the 

300 period from January 1979 to December 2007, inclusive; no corresponding data were available 

301 prior to November 1978. Both ice variables to some degree measure the amount of open water 

302 over the shelf, one directly (area) and one indirectly (extent). Ice extent in the Ross Sea region 

303 varies with winter/spring wind strength and persistence, as does the extent of the coastal Ross 

304 Sea polynya (explained more fully in Ainley et al. 2005, 2010): more wind pushes the ice farther 
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305 north, but also frees the coastal area of ice. It there is little wind, there is little water among ice 

306 floes and ice area is maximized. 

307 Statistical analyses were primarily conducted on average monthly values of ice area and 

308 ice extent. The maximum and minimum values for the 12 monthly averages were determined 

309 for each year. The minimum was always in February, but the maximum could occur in August, 

310 September, October, or November. In addition, we determined daily maximum and minimum 

311 values. However, results were not improved using daily rather than monthly minimum or 

312 maximum values, and hence we show only results using the month-specific values, whether by 

313 specified calendar month, e.g., "September", or with reference to maximum or minimum 

314 monthly value. 

315 Analyses were also conducted on duration of the "ice season" or "open water season": 

316 (1) duration (number of days) that ice area exceeded an upper threshold value (3,037,000 km 2 ) 

317 relative to the total number of days for which the metric was calculated for that year, and (2) 

318 duration (number of days) that ice area was below a lower threshold value (628,000 km 2 ) 

319 relative to the total number of days calculated for that year. The upper threshold was chosen to 

320 be the 95 th percentile value of ice area in the year with the lowest maximum value (1986). The 

321 lower threshold was chosen to be the 5 th percentile value of ice area in the year with the 

322 highest minimum value (2003). However, these measures of duration above an upper threshold 

323 or below a lower threshold were no better predictor variables with respect to toothfish metrics 

324 than were the maximum or minimum monthly averages of ice area or ice extent explained 

325 above. Thus, the results with respect to relative ice season duration are not shown. 
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326 Additional ice metrics analyzed included: timing (day of year) of minimum and 

327 maximum, the difference in ice area or ice extent between minimum and the following 

328 maximum, difference between the minimum and the previous maximum in terms of ice area or 

329 ice extent, the difference in days between the minimum and the maximum, and the rate of 

330 change (difference in ice area or extent between maximum and minimum divided by number of 

331 days between minimum and maximum). Again, none of these variables proved superior in 

332 predicting fish metrics and so results with respect to these variables also are not shown. 

333 

334 Other environmental variables. We examined the importance of several environmental 

335 variables that by influencing ocean characteristics may in turn directly or indirectly affect 

336 toothfish, as they have other Ross Sea mesopredators (Testa et at. 1992, Wilson et at. 2001, 

337 Ainley et al. 2005, Rotella et al. 2009). These other environmental variables were the Southern 

338 Oscillation Index (SOI), the Southern Annular Mode (SAM; also called Antarctic Oscillation 

339 Index), and a measure of the size of the Ross Sea Polynya. It is well known that the state of the 

340 sea ice in this part of the Southern Ocean is affected by the atmospheric processes embedded 

341 in SOI and SAM (Jacobs 2006, Stammerjohn et al. 2008). 

342 We used monthly SOI to calculate six metrics. The first three were: average annual SOI 

343 (January to December), the average for January to June, and the average for July to December. 

344 The other three were quadratic transformations of the corresponding SOI metric, designed to 

345 measure increasing effects of SOI associated with strong El Nino events, during which the SOI is 

346 strongly negative. The transformations were of the form (SOI- c) 2 where c = maximum (positive 
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347 value) of the SOI metric (12-mo average or either of the 6-mo averages). (SOI - c) was thus 

348 always zero or negative and the stronger the El Nino the more negative was the value. Squaring 

349 (SOI - c) resulted in all non-negative values, with strong El Nino being associated with the 

350 largest positive values of the transformed metric and the strongest La Nina in the time series 

351 getting a score of zero (following the approach of Lee et al. 2007). In no case was a second- 

352 order polynomial (i.e., with 2 df) superior in terms of AIC to the respective quadratic 

353 transformation of SOI (i.e., with 1 df) when analyzing fish condition or length. 

354 For all single-variable regression models, we report the regression coefficient (b, which 

355 is also commonly symbolized as p). Where more than one environmental variable was analyzed 

356 in the same model, we report the multiple regression coefficient (b) for the effect of a variable 

357 while controlling for all other variables in the model, symbolized. To aid in the comparison 

358 among variables (measured on very different scales, e.g., SOI and maximum ice area), we also 

359 present the standardized multiple regression coefficient (which we refer to as “beta"), which 

360 like a correlation coefficient, scales from -1 to +1 (Kutner et at. 2005). 

361 SAM was represented by the value on 1 January each year (from 

362 http://www.antarctica.ac.uk/met/gjma/sam.html; Fig. 2) and is a proxy for strength of the 

363 Southern Ocean westerlies, which in turn affect the formation, extent and persistence of sea ice 

364 (Marshall et al. 2003, Russell et al. 2006, Stammerjohn et at. 2008). Using SAM values from 1 

365 January, about 9 to 11 months before toothfish were captured as part of this study, means that 

366 a lag of -1 was used; lag of 0 means SAM is calculated in the same calendar year as applies to 
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367 the sample of toothfish; lag of 1 means SAM is calculated between 1-plus years (21 to 23 

368 months) before the toothfish sample was obtained, etc. 

369 Among these environmental variables, we examined condition, length, and numbers 

370 caught per unit effort in relation to same-year values, as well as with lags of 1, 2, 3, 4, and 5 yr. 

371 One rationale for examining lags was that condition and numbers of toothfish could well reflect 

372 availability of their prey, i.e., Antarctic silverfish. For example, in McMurdo Sound toothfish 

373 consume silverfish that are 81-195 mm SL with modal values of 100-150 mm (Eastman and 

374 DeVries, unpublished data) and, on the basis of Hubold and Tomo's (1989) aging studies for 

375 Weddell Sea silverfish, these specimens are 4+ years old. Thus, we reasoned ice (or other 

376 environmental) conditions in a given year may influence availability of prey 1 to 5 years later. 

377 

378 

379 Results 

380 Characterization of the catch 

381 Between 1972 and 2010, 5587 toothfish were caught, including 835 kept for experiments 

382 and/or whose survival was in doubt. The largest fish caught was 193 cm (92 kg) and the 

383 smallest was 81 cm (4.5 kg); individuals between 120 and 170 cm predominated, with very few 

384 <90-100 cm (Fig. 3). 

385 Of the 4752 fish tagged and released in the study, 17 have been recaptured (0.4%) with 

386 the annual growth rate being 2.0 cm in TL and 1 kg in mass per year (DeVries, unpubl. data). 
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387 These growth rates are slightly below those reported from analyses of the industrial catch (cf. 

388 Horn 2002, Brooks et al. 2010). Most of the recaptures (12) occurred at the McMurdo fishing 

389 site, but one was recaptured by the Russian vessel, Ycmtar, and four have been reported by the 

390 New Zealand fleet (DeVries, unpubl. data). One tagged individual was recaptured >1300 km 

391 north of McMurdo Sound, indicating migration north of the Ross Sea, possibly for spawning 

392 (Hanchet et al. 2008, 2010). Most recaptured McMurdo fish were encountered again 4-5 yr 

393 after release, with the longest interval being 18 yr. The recapture rate in the fishery during the 

394 most recent years, 2004-2009, for a select number of boats and not including fish recaptured 

395 the same year, is 4.96% (CCAMLR 2010). The discrepancy between our low recapture rate and 

396 this one could well be due to 1) McMurdo Sound fish being more transient than those on the 

397 main fishing grounds (continental slope) and/or 2) the fact that recaptures of our fish by the 

398 fishery were under-reported, with many vessels other than those from NZ (just 4 of 15-20 

399 permitted) reporting no recaptures at all (CCAMLR 2006). This would especially be true in early 

400 years of fishery, when no industry tagging program was in operation; and of course no data 

401 come from the Illegal, Unreported and Unregistered fleet. One of our tags was found on a 

402 vessel's processing room floor, indicating lack of vigilance. 

403 

404 Fish Length 

405 Annual variation in total length, considered an approximate proxy for age, was analyzed with 

406 respect to three metrics: 50 th , 75 th (i.e., upper quartile), and 95 th percentiles. In general, the 

407 overall pattern during the study period was similar for each of these metrics. 
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408 The most parsimonious model, favored by AIC, describing the pattern of change in 

409 length across years was a two-part spline with change point (i.e., "knot") at 1992, with 

410 increasing TL to 1992 and decreasing TL subsequently (Table 2, Fig. 4). These models were 

411 superior to a single linear trend, quadratic polynomial, or cubic polynomial model, or other 

412 simple transformation, as determined by AIC, and they accounted well for the pattern 

413 observed, except for 1972 being an outlier with respect to the predicted pattern. In regard to 

414 the 75 th percentile, for the linear trend 1972-1992, b = 0.580 (± 0.123 SE, t = 4.7, P<0.0001) and 

415 for the linear trend 1992-2010, b = -0.526 (± 0.242, t = -2.2, P<0.04). The change in trend was 

416 significant (Table 2). In regard to the 95 th percentile, for the line segment 1972-1992, b = 0.534 

417 (± 0.169, t = 3.6, P = 0.004) and for the line segment 1992-2010, b = -0.662 (± 0.331, t = -2.0, P = 

418 0.057). The change in slope for the 95 th percentile was also significant (Table 2). Finally, the 50 th 

419 percentile (not shown) showed the same pattern, and was statistically significant with respect 

420 to analogous measures. 

421 In order to assess whether or not the length statistics of the small sample of fish caught 

422 during 1991 and 1991 were mainly responsible for the observed 2-part trends, we also analyzed 

423 trends using individuals rather than years, analyzing with and without data from 1991 and 

424 1992. Models included effect of calendar month (analyzed as categorical variable; Table 3). The 

425 2-part spline model including 1991 and 1992, as expected, was highly significant (Table 3), as 

426 was the model excluding 1991 and 1992 data [P <0.0001, Table 3). The coefficient, b, for the 

427 1972-1992 slope was 0.362 (SE = 0.041), and that for the 1972-1990 slope was 0.351 (SE = 

428 0.041; in both cases, P <0.0001). Similarly, the coefficient, b, for the 1992-2010 slope was - 

429 0.291 (SE = 0.104) and that for the 1997-2010 slope was -0.272 (SE = 0.104; P = 0.005 and P = 
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430 009, respectively). Finally, the change in slope from the earlier period (up to 1992) to the later 

431 period (from 1992) was significant whether or not 1991 and 1992 were included (Table 3). 

432 

433 Fish Condition 

434 Condition (K) of the Antarctic toothfish caught varied by month. K was relatively high in October 

435 and November, but lower in September and December (Fig. 5). Variation among months was 

436 significant (F 3)5 394 = 4.09; P = 0.007); the pattern and significance of monthly variation in 

437 condition were very similar after controlling for annual variation (year as factor). 

438 Across all years, there was no significant linear trend in toothfish condition (Fig. 6). 

439 Instead there was a trend for mean K to increase from 1972 to 1991 and/or 1992 and decrease 

440 subsequently. Mean K was higher for 1991 and 1992 than in any other year; nevertheless, an 

441 increasing trend is evident from 1972 to 1990. The pattern evident indicated that change point 

442 analysis, with either 1991 or 1992 being the knot, would be most appropriate. Adjusted R 2 was 

443 maximum and AIC was optimized when the change point was at 1992 rather than at 1991 

444 (Table 2). For the line segment 1972-1992, b = 0.00358 (± 0.00094, t = 3.8, P = 0.001), and for 

445 the line segment 1992-2010, b = -0.00927 (± 0.00180, t = -5.2, P <0.001); the difference in 

446 slopes of the two lines was significant (Table 2). 

447 AIC for a quadratic model for year was poorer than for the change-point analysis. More 

448 specifically, there is no evidence that the trend from 1972 to 1990 was anything but linear 

449 (linear trend through 1990: P <0.0001, quadratic term for that time period, P >0.2) and. 
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450 therefore, the change-point analysis is preferred, with a change in trend either in 1991 or in 

451 1992 (Figure 6). 

452 In order to assess whether or not mean condition of the small sample of fish caught 

453 during 1991 and 1991 were mainly responsible for the observed trends in condition, we also 

454 analyzed trends using individuals rather than mean values, analyzing with and without data 

455 from 1991 and 1992. As with length, models included effect of calendar month (analyzed as 

456 categorical variable; Table 3). The 2-part spline model including 1991 and 1992, as expected, 

457 was highly significant (Table 3), as was the model excluding 1991 and 1992 data (Table 3). The 

458 coefficient, b, for the 1972-1992 slope was 0.00338 (S.E. = 0.00039), and that for the 1972-1990 

459 slope was 0.00319 (S.E. = 0.00039 ;in both cases, P <0.0001). Similarly, the coefficient, b, for the 

460 1992-2010 slope was -0.00938 (S.E. = 0.00099 ) and that for the 1997-2010 slope was -0.00910 

461 (S.E. = 0.00099 ;P <0.001 in both cases). Finally, the change in slope from the earlier period (up 

462 to 1992) to the later period (from 1992) was significant whether or not 1991 and 1992 were 

463 included (Table 3). 

464 

465 Axe-handle fish. The proportion of axe-handle fish decreased during the period 1972 to 

466 1992 from about 3 - 7% in the mid-1970's to 0 - 3% in the period 1987 to 1992, and then 

467 increased from 1992 on (Fig. 7). The one exception was 2010, in which there were no axe- 

468 handle fish, but only 7 total fish were caught that year. Among the set of models examined, the 

469 preferred model (as determined by AIC) was a two-part spline, with linear decrease up to 1992 

470 (b = -0.00167 ± 0.00076, t = -2.19, P = 0.039), and a linear increase since (b = 0.00420 ± 0.00145, 
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471 t = 2.90, P = 0.008). The linear trend in each time period was significant (see below) and the 

472 trends differ significantly from each other (Fi /2 4 = 8.68, P = 0.007). 

473 

474 Catch per Unit Effort 

475 Catch rate varied by month within each season (Fig. 5). The lowest capture rate occurred in 

476 September (even lower in August but only two days of fishing were available), but it then 

477 increased dramatically in October. Rates for November and December remained high, though 

478 with a slight downward trend from November to December. There was no significant variation 

479 in catch rates among October, November, December (P > 0.4); but each of these differed from 

480 the September catch rate (P < 0.001). 

481 In regard to annual trends, unlike condition and length, there was no peak in CPUE 

482 midway during the study period. Instead, captures were fairly stable through 2001 (Fig. 8), with 

483 a slight, non-significant decline in capture rate evident late in this period, equivalent to a 

484 decline in CPUE of 1.24% per year (back-transformed geometric mean SE = 1.76%: analysis of In- 

485 transformed monthly capture rates, controlling for month of capture, slope not different from 

486 zero, P > 0.4). The pattern then changed dramatically, but unfortunately only 3 years of effort 

487 data were available after 2001. For those 3 years, captures per 10 hooks deployed, summarized 

488 per year, ranged 0 - 0.4 fish, compared to 1975-2001 when the range was 1.0 - 3.4 fish. Though 

489 the information is not included in our analyses, fishing was undertaken by another researcher 

490 (G. Flofmann) in the same localities in 2006 for 4 weeks: no fish were caught but effort data are 


491 


not available. 
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492 Realizing that data were available for only 3 of the 9 years after 2001, and that these 3 

493 years occurred at the end of this latter period, the best-fitting, most parsimonious model 

494 (optimizing AIC) indicated a step change to a different (and much lower) level of capture rate. 

495 This, however, must be viewed with caution with so few years since 2001. The pattern depicted 

496 is preferred by AIC when compared to the full set of competing models: (1) no change, (2) two 

497 linear segments joined at 1992, (3) two linear segments joined at 2001, (4) two linear segments 

498 joined at any year between 1992 and 2001, (5) two parallel lines with a step change occurring 

499 at or after 2001 (i.e., slope is non-zero but the same across all years; only the elevation changes 

500 at or after 2001), (6) inverse transformation of year, (7) year as a quadratic function, and (8) 

501 year as a cubic function. Nevertheless, with only three years of catch per unit effort available 

502 post-2001, it is not possible to adequately estimate the post-2001 trend. Essentially, all that can 

503 be said is that fishing success became dramatically reduced after 2001. 

504 

505 Influence of environmental variables on fish metrics 

506 Fish Length. Exploratory analyses of the effects of the full set of candidate ice variables 

507 (see Methods) on our three TL metrics revealed four of particular interest: minimum ice area, 

508 minimum ice extent, maximum ice area, and maximum ice extent. Each of these was analyzed 

509 with respect to lags of 0 - 5 yr, yielding 24 ice variables (= 4 variables x 6 lags). For 95 th 

510 percentile length, none of the 24 ice variables was significant. For 75 th percentile length two 

511 variables were significant: maximum ice area with 4-yr lag and maximum ice area with 5-yr lag 

512 (Fig. 9); for 4-yr-lag, weighted regression, b = -0.0000084 (± 0.0000039, t = -2.15, P = 0.048; 
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513 Table 4). We do not show results for 5-yr-lag because this relationship, though significant (b = - 

514 0.0000089 ± 0.0000041, P = 0.049), is driven entirely by results from 1991 and 1992. 

515 Comparison of correlation coefficients between 75 th percentile length and each of the 

516 two ice variables (maximum ice area with 4-yr lag vs 5-yr lag) is instructive, with and without 

517 1991 and 1992. For the 5-yr lag, the correlation coefficient with respect to 75 th percentile 

518 length was r =-0.499, including the 2 years, but r = -0.050 without those 2 years. For the 4-yr 

519 lag, the correlation coefficient with respect to 75 th percentile length was r = -0.486 including the 

520 2 years, but only dropped to r = -0.404, excluding the two years. Thus, the relationship 

521 between 75 th percentile length and maximum ice area, with 4-yr-lag, does not hinge on these 2 

522 years having unusually large upper quartile values. 

523 Median length (50% percentile) also demonstrated a significant relationship to 

524 maximum ice area 4 years previously (Table 4; Fig. 9), and one that was stronger than that 

525 shown for upper quartile length (for median length: correlation coefficient, r = -0.608, P < 0.01; 

526 compared to r = -0.486 for upper quartile length and the same ice variable). This relationship is 

527 robust with respect to values for any individual year. Whereas 1991 is an influential 

528 observation, 1992 in this case is not (median length in 1992 = 136, compared to median length 

529 for entire sample = 134). Excluding 1991 still yields a relationship that is close to significant (P = 

530 0.075). 

531 

532 Fish Condition. Two distinct ice variables provided the best statistical predictors of 

533 changes in K\ minimum ice area (always in February) of the current year (8-10 mo prior to 
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534 fishing season) and maximum ice extent 3 yr previously (Fig. 10). While minimum ice area could 

535 also be a proxy for amount of open water, understanding maximum ice extent is far more 

536 complex, owing to the fact that large scale ice extent is a function largely of wind strength and 

537 persistence in the Ross Sea sector. Stronger winds (greater ice extent) also usually mean that 

538 the coastal polynyas are more prevalent (though not always, depending on temperature), a 

539 subject we shall return to in the Discussion. Regardless, these two variables (minimum ice area, 

540 maximum ice extent 3 yr previous) were relatively uncorrelated with each other (r = 0.324). 

541 There were several other variables (such as length of ice season) that provided good predictive 

542 ability for explaining change in condition, though not quite as high, but also highly correlated 

543 with the two primary variables. Thus, as noted in Methods, we keep those results in mind but 

544 do not report them. 

545 The relationship between toothfish condition and minimum ice area 8 - 10 mo earlier 

546 was non-linear: mean K declined with increasing minimum ice area (less open water than 

547 usual), but in an accelerating (down-turned) fashion. In a weighted regression, mean K in 

548 relation to minimum ice area, the relationship best reflected a quadratic transformation of 

549 minimum ice area (predictor variable = (minimum ice - c) 2 where c = smallest minimum ice area 

550 observed in any year = 131,448 in 1979) and was significant: b = -2.39 x 10 13 (± 1.07 x 10 13 , t = - 

551 2.23, P = 0.039; Table 4). 

552 When analyzing the data at the individual level, the effect of minimum ice area 

553 (quadratic transformed) on condition, while controlling for calendar month of capture, was also 

554 significant, b = -1.52 x 10 13 (± 4.08 x 10 14 t =-3.73, P = 0.0002). Furthermore, even excluding 
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555 the influential year 2003 (greatest minimum ice area and lowest mean condition index in the 

556 data set), the result of analyzing individual fish condition demonstrated a significant negative 

557 effect of minimum ice area (quadratic transformed), b = -1.14 x 10 13 (± 4.24 x 10 " 14 , t =-2.69, P = 

558 0.007). Minimum ice area was a better predictor (as determined by AIC) of mean condition than 

559 was date of minimum or minimum ice extent. 

560 The second primary ice variable to account for annual variation in mean condition was 

561 maximum ice extent 3 years earlier (Fig. 10); other lags showed no significant relationship, (0, 1, 

562 2, 4, or 5). Maximum ice extent was a better predictor than was ice extent in any specific 

563 month (e.g., October). In this case the relationship was linear. In a weighted regression, mean K 

564 in relation to maximum ice extent with 3-yr-lag was: b = -7.18 x 10 " 8 (± 3.19 x 10" 8 , t = -2.23, P = 

565 0.04; Table 3). 

566 There was no significant relationship (P >0.15) between proportion axe-handle fish and 

567 minimum ice area in the same year nor between proportion axe-handle fish and maximum ice 

568 extent 3 years previously, the two variables that correlated strongly with mean condition. 

569 However, there was a significant relationship between the proportion of axe-handle fish and 

570 maximum ice area 4 years previously, the same variable that predicted 50 th percentile and 75 th 

571 percentile length (Fig. 7). With increasing maximum ice area, the proportion of axe-handle fish 

572 increased 4 years later (adj r 2 = 0.2698, F 145 = 6.91, P = 0.019). 

573 

574 Change in Catch Rate (CPUE). There was no evidence that any environmental variable 

575 (ice variables, whether lagged or not, or SOI variables) affected toothfish CPUE (In-transformed) 
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576 or annual changes in toothfish CPUE. Furthermore, there is no significant or AlC-justified 

577 association between CPUE (or changes in CPUE) and any of the ice variables, lagged or not, 

578 during the pre-crash period, 1972 to 2001. As noted, despite environmental variability, there 

579 was little change in capture rates from one year to the next up to 2001 (Figure 8). 

580 

581 Relationships to the Southern Oscillation. No significant relationship existed between 

582 SOI and the three types of fish metrics analyzed. This is not to say that SOI is unimportant, only 

583 that the effect of SOI is not direct. Instead, SOI during the months July-December is strongly 

584 associated with maximum ice area in the same year (this reaffirms previous analyses, i.e., 

585 Wilson et al. 2001); b = 108,700 (± 2400, t = 4.52, P = 0.0001; Fig. 11). A linear relationship 

586 between the two variables is best supported and confirmed by AIC (Table 3). 

587 Because maximum ice area reflected SOI in July to December and maximum ice area 

588 predicted 75 th percentile length of toothfish 4 years later (see above), there was consequently a 

589 modest, non-significant, negative correlation between SOI in July-December and 75 th percentile 

590 length 4 years later (r = -0.255, P > 0.2; Table 4). Flowever, once we accounted for the effect of 

591 maximum ice area on length (with 4-yr-lag), the effect of SOI in July-December on 75 th 

592 percentile length 4 years later was even less significant and reversed sign (partial regression 

593 coefficient, b = 0.585 ± 0.721, t = 0.81, P > 0.4; standardized partial regression coefficient beta = 

594 0.225.), while the effect of maximum ice area on length 4 years later remained, even while 

595 controlling for SOI in July-December (partial regression coefficient, b = -0.0000107 ± 

596 0.00000485, t = -2.21, P = 0.044; standardized partial regression coefficient beta = -0.613). 
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Results were very similar to that found for 75 th percentile length analyzing either 50 th or 
95 th percentile length in relation to the same two independent variables: SOI (July-Dee) with a 
4-year lag and maximum ice area with a 4-year lag. The effect of maximum ice area was 
significant, controlling for SOI, while the effect of SOI was not. 

It was also the case that SOI influenced maximum ice extent and thus exerted an 
indirect effect on mean condition. However, just as with 75 th percentile length, there was no 
direct effect of SOI on condition 3 years later, once we controlled for maximum ice extent 
(partial regression coefficient b = 0.0063 ± 0.0082, t = 0.73, P > 0.4, standardized partial 
regression coefficient beta = +0.217). Note that SOI in the current year was not correlated with 
mean condition or any of the length metrics, nor was SOI associated with minimum ice area in 
the current year. 


Relationships to Southern Annular Mode. There was little evidence of a direct effect of 
SAM (measured on a continuous scale) on toothfish length or condition, whether lagged or not. 
The regression of mean condition or the three percentiles of length on SAM was not significant 
in every case, with lags of 1 to 5 years. 

However, the sign of SAM changed from predominantly negative prior to 1976, to 
increasingly positive from 1976 to the early 1990s, and then to mostly positive thereafter; 
during the period from 1972 (first year of this study) to 1992, 12 out of 21 years were < 0, 
whereas from 1993 to 2011, 15 out of 19 years were positive (LRS = 5.60, df = 1, P = 0.018; Fig. 
2). In accord with the increasing windiness as a result of positive SAM, sea ice extent (and 
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618 related aspects of sea ice) changed as well (Fig. 11). During the more "negative" period of SAM 

619 (1972 to 1992), condition and length both increased. During the (mostly) "positive" period of 

620 SAM (1993 to 2011), condition and length both declined significantly. 

621 While there was no relationship between condition or the length metrics and SAM as a 

622 quantitative variable, 75 th percentile length was significantly greater in years in which SAM was 

623 positive 4 years earlier as compared with years in which SAM was negative (<0) 4 years earlier: 

624 145.85 (± 1.17; n = 14 years) vs. 142.14 (± 1.20; n = 13 years; P = 0.036; results of weighted 

625 regression). Results were even stronger comparing 50 th percentile lengths. The median fish 

626 length 4 years after a year with positive SAM was 136.39 (SE = 1.21) compared to median fish 

627 length of 131.84 (SE = 1.13) 4 years after a year with SAM of 0 or less ( P = 0.011). 

628 

629 Discussion 

630 In this study we found that (1) the length-frequency of the McMurdo Sound catch was 

631 dominated by fish >100 cm and in the upper two-thirds of the overall distribution exhibited in 

632 the industrial toothfish catch; (2) McMurdo Sound fish length and condition increased from the 

633 early 1970's to the early 1990s and then decreased, the proportion of "axe handle" fish (poor 

634 condition) in the catch changing inversely with mean fish condition; (3) fish length and 

635 condition decreased with increasing sea ice cover and/or inversely with the extent of open 

636 water, lagged by 8 mo to 5 yr; and (4) catch per unit effort, controlling for within-season 

637 variability, was constant through 2001 and then decreased dramatically. Below, we propose 

638 and discuss hypotheses explaining these patterns. 
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639 

640 Toothfish length-frequency 

641 Models to estimate stock size based entirely on what can be gleaned from the industrial catch, 

642 as is the case for the Ross Sea region toothfish fishery (CCAMLR 2005-2010), assume that after 

643 making certain statistical assumptions and corrections, e.g., by virtue of mark-recapture data, 

644 models are sufficient to inform the management strategy. One main assumption is that the 

645 size-age distribution in the catch, with adjustment, is representative of the full population. The 

646 disparity in the size distribution in the McMurdo Sound catch compared to the fishery catch is 

647 therefore of concern. The McMurdo Sound catch, ranging from 81-193 cm with a broad mode 

648 at 120-170 cm, reasonably mirrors the upper two-thirds of the commercial catch, i.e. "...50 to 

649 180 cm... .In all seasons, there was a broad mode of adult fish at about 120-170 cm..." (CCAMLR 

650 2010, p. 6); ~15-20% of the industrial catch overall appears to be offish <100 cm, and in years 

651 when there are more of these smaller fish the conclusion is that the fishery centered more over 

652 the Ross Sea shelf than in other years (CCAMLR 2010). In fact, with very few fish <100 cm (2.5%) 

653 in the McMurdo Sound catch, this disparity is even more apparent when compared to the size- 

654 frequency distributions in Hanchet et al. (2010). In the latter, only fish captured by the fishery 

655 over the shelf are displayed and >50% are <100 cm, with the fishery expending significant effort 

656 over the shelf (though less than on the slope) including waters close to the scientific fishing 

657 location (i.e., on the opposite side of Ross Island (Hanchet et al. 2008). 

658 There are at least two factors that could explain the discrepancy. First, owing to lack of 

659 buoyancy in fish <100 cm (Eastman and Sidell 2002, Near et al. 2003), the vertical setline failed 
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660 to attract many of the smaller fish that remain on the bottom (Eastman and Barry 2002, Near et 

661 al. 2003), whereas the industrial longlines, deployed along the bottom, caught these fish in 

662 greater proportion. Second, the discrepancy could be due to the fact that large fish ate most of 

663 the smaller fish caught on the vertical set line. However, of the stomachs inspected in the 

664 McMurdo catch [n = 58), evidence for cannibalism was absent (Eastman 1985a, b; La Mesa et 

665 at. 2004). That large toothfish, at least under ice cover, reside high in the water column, as 

666 indicated by the catch we report herein, is a fact confirmed by crittercams placed on Weddell 

667 seals as well as on a number of occasions by ROVs (to within 12 m of the surface: Fuiman et al. 

668 2002, Kim et al. 2011); that small toothfish do occur in McMurdo Sound is indicated by the 

669 catch of Weddell seals, who take the entire size range (cf. Ponganis and Stockard 2007, Ainley 

670 and Siniff 2009). In any case, it would appear that assumptions made in CCAMLR's fishery 

671 models about the representativeness of age-size structure as revealed by longlines need to be 

672 reconsidered. 

673 In the context of the discrepancy between vertical and horizontal-bottom sampling it is 

674 interesting to consider why this species is one of the few neutrally buoyant Southern Ocean 

675 fish. The most parsimonious explanation is probably historical. Given the eradication of the 

676 previous taxonomically diverse Eocene fauna, the modern notothenioids diversified into a 

677 developing ecosystem with vacant water-column niches. As inferred from molecular sequence 

678 data, these niches were filled first by the neutrally buoyant clade about 24 million years ago 

679 (Near 2004). The water column contained few fishes and underexploited resources were 

680 available— the Antarctic toothfish diversified coincident with silverfish in response to this 

681 unusual situation. The questions remain: Is neutral buoyancy the means to exploit the high 
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682 prevalence of silverfish (another neutrally buoyant species), the means to escape cannibalism, 

683 or both (Pinkerton et al. 2007)? Is the lack of any water-column species in the diet of longline- 

684 caught fish (Fenaughty et al. 2003) a reflection of benthic longlines catching mainly fish that 

685 reside on or near the bottom? Why, supposedly, are toothfish only found mostly on the bottom 

686 in ice-free waters? Regardless, it appears that the clear ontogenetic shift in Antarctic toothfish 

687 from dwelling on the bottom as small fish to somewhere in the water column as large fish, at 

688 least in ice covered waters, is a major aspect of this species' natural history (DeVries and 

689 Eastman 1981, Near et ol. 2003), with further implication of whether or not removing large fish 

690 results in the increase of small fish (through competitive release), i.e. the basis for modern 

691 fishery theory (cf. Constable et al. 2000, Longhurst 2010). In the hypothetical life history 

692 scenario proposed by Hanchet et al. (2008) no mention is made of the ontogenetic shift in 

693 buoyancy and its implications, other than to say that as fish mature they occupy waters having 

694 deeper bathymetry. The assumption is made that this species is entirely bottom dwelling 

695 (except perhaps larvae). 

696 One further interesting question that arises from the size-frequency of McMurdo Sound 

697 fish compared to that of the industrial catch revolves around the age-related fecundity of this 

698 species. All large McMurdo Sound fish are sexually mature, but in gonadal resting stage; some 

699 females showed signs of previous spawning (Eastman and DeVries 2000). If the industrial catch 

700 underestimates the prevalence of large fish over the shelf, it would seem that estimates of 

701 age/size at maturity, which is based on the size-frequency seen in the industrial catch, needs to 

702 be revised once more (recently it was doubled from ~8 to ~16 y; Parker and Grimes 2010). 
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703 

704 Trends in toothfish length and condition 

705 Both total length and condition of toothfish in the McMurdo Sound catch increased from the 

706 mid 1970s to the early 1990s, and then switched to gradually decrease. Because of the low 

707 sample size for 1991 and 1992, and the lack of any data during 1993-1996, it is difficult to say 

708 when length and condition trends actually leveled off in the early to mid-1990's before 

709 decreasing. However, in the period up to and including 1990 there was no indication of leveling 

710 off. Instead, a linear trend was evident in condition and length from 1972 to 1990. Whether or 

711 not the decreasing trend is represented in the industrial catch, which only began (legally) in 

712 1996 (CCAMLR refers this as the 1997 fishing season), we do not know. Seemingly this would be 

713 difficult to ascertain in a sampling regime that is not only non-random but is geared to catch the 

714 largest fish as fast as possible before the season closes (Constable et at. 2000, Brooks 2008). 

715 Regardless, it appears that growth rate in this species is not a constant but is subjected 

716 to environmental and potentially anthropogenic factors. According to fishery theory, as fish 

717 (especially large ones) are removed from the stock, fish size and condition should show a 

718 positive growth as smaller fish are released from competition. This is rarely demonstrated 

719 (Longhurst 2010) and is not evident in the McMurdo Sound data. What then are the factors 

720 that could account for the observed increases and then decreases in length and condition? If 

721 change in length is a function both of food availability and age, was there 1) a release from 

722 predation during the 1970-80s, thus allowing older fish to remain in the population longer, 

723 and/or 2) more favorable foraging related to competitive release or climatic factors that 


724 


influenced silverfish? 
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725 

726 Release from predation. The southern Ross Sea was the scene for a number of "heroic" 

727 expeditions during the early half of the 1900s, and they all took large numbers of Weddell seals 

728 for human and dog food. The southern Ross Sea seal population was thought to have recovered 

729 by the 1960s (Stirling 1971), but then the NZ Antarctic Programme began to kill 50-100 seals 

730 annually to feed sled dogs and by the time this ended in the mid-1980s, the population, at least 

731 as measured in McMurdo Sound was halved, from 3000 to 1500 seals. Then over the next 

732 decade, it slowly recovered (partially) to ~2000+, where it has remained ever since 

733 (summarized in Ainley 2010). Given that Weddell seals eat a significant number of Antarctic 

734 toothfish (Ainley and Siniff 2009, and references therein; Kim et al. 2011), this marked decrease 

735 in Weddell seals in the southern Ross Sea is consistent with a lessening of predation pressure 

736 on Antarctic toothfish. Certainly, the lag that would occur between these two long-lived, late 

737 maturing species, predator versus prey, would be large. Added to this is the fact that another 

738 toothfish predator, sperm whales (admittedly at the periphery of their range), owing mainly to 

739 contraction of range as whalers depleted the stocks in the warmer parts of the Pacific during 

740 1700-mid 1900s, even now are a shadow of their former numbers in the Ross Sea sector 

741 (Whitehead 2000). Their numbers in the Ross Sea sector also may well have begun to recover in 

742 recent decades as they have elsewhere (see for instance, Whitehead et al. 1997, Branch and 

743 Butterworth 2001). The same is true for southern elephant seals, still another toothfish 

744 predator, which disappeared from the Ross Sea during the 1980s, their foraging range 

745 contracting as the Macquarie Island breeding population became severely reduced 

746 (summarized in Ainley 2010). 
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747 

748 More favorable foraging related to competitive release or climatic factors. The industrial 

749 take of minke whales during the 1970s-early 1980s (~20,000 taken from Ross Sea region), which 

750 appears to have resulted in the population increase of trophically competing Adelie Penguins at 

751 the time (Ainley et at. 2007), also could have benefited Antarctic toothfish, both the whales 

752 (and penguins) and the toothfish being significant predators of Antarctic silverfish (Ballard et at. 

753 2011). Furthermore, the time of the switch in toothfish length and condition trends from 

754 positive to negative corresponds to when the minke whale population was deemed to have 

755 reached recovery (roughly 1986-1990; Branch 2006, Ainley 2010). 

756 There is also the possibility that the shifting climate of the Ross Sea region has affected 

757 (is affecting) the abundance of Antarctic silverfish, a subject about which we know even less. 

758 According to La Mesa et al. (2010), the prevalence of larval silverfish in the western Ross Sea 

759 was an order of magnitude higher in the 1999-2000 austral summer, when the Ross Sea 

760 polynyas opened early and large, compared to the austral summers of 1997-98 and 2003-04; 

761 with more open water, earlier, primary production is increased and, hypothetically, silverfish 

762 survival is enhanced owing to increased food. Silverfish spawn during the winter and their eggs 

763 can be found among the frazil ice associated with heavy pack or fast ice (Vacchi et al. 2004). 

764 Perhaps this is a protective strategy, and more ice at this time is beneficial, but once the eggs 

765 hatch the larvae need to forage on microbial organisms. In any case, though most of the 

766 silverfish contained in toothfish stomachs in McMurdo Sound during late spring are 100-150 

767 mm SL (Eastman and DeVries, unpublished data) and on the order of 4+ years old (Hubold and 

768 Tomo 1989), we did find that toothfish condition is positively affected by extensive open water. 
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769 and hence water open to insolation, during mid-summer. We also found the effect of a 3-yr-lag 

770 in ice extent on toothfish condition, and note that Ross Sea primary productivity cycles on a 2-4 

771 year period owing to a combination of nutrient and sea ice factors (Peloquin and Smith 2007). 

772 Consistent with the 4 yr and older age classes of silverfish that toothfish eat in mid-summer, we 

773 also found a negative relationship between toothfish length and maximum sea ice area 4-5 

774 years earlier. Greater ice area means less open water. That's not to say that we understand the 

775 interconnections, but these few factors and relationships would seem to implicate the 

776 availability (and possibly condition) of silverfish as being important to toothfish growth and 

777 condition. Obviously, we need to learn considerably more, not just about toothfish, but about 

778 silverfish as well, if we are to understand toothfish. Making what we can of the trends in axe 

779 handle fish, we note the consistency between the decreasing trend in prevalence while 

780 toothfish length and mean condition were increasing during the first half of the study period, as 

781 opposed to the increased prevalence of axe handles as length and condition began to decrease. 

782 Another factor to consider is the Southern Annular Mode, a major driver of climate in 

783 the Southern Ocean (Thompson and Solomon 2002, Stammerjohn et at. 2008). During the 

784 period when toothfish length and condition were increasing, beginning in the mid 1970s SAM 

785 was increasingly becoming less negative and more positive until becoming almost entirely 

786 positive starting at the inflection point of both toothfish length and condition. Therefore, owing 

787 to increased windiness blowing south to north (as well as westerly winds in the north increasing 

788 Ekman transport), sea ice extent was expanding, early appearance of coastal polynyas was 

789 becoming more reliable, and the sea ice season was lengthening (Parkinson 2002, Stammerjohn 

790 et al. 2008). It would seem that more reliably present coastal polynyas in the Ross Sea would 



38 


791 favor silverfish production, but the decreasing annual period of open water would not. Other 

792 than the correspondence in the inflection points of length and condition and the onset of 

793 continually positive SAM, we again emphasize that we don't understand much about this 

794 predator-prey system, toothfish and silverfish. The correspondence in trends between those 

795 changes in trend and SAM, though, is noteworthy and deserving of closer study. 

796 That the Southern Oscillation Index had no direct correlation to variation in any 

797 Antarctic toothfish parameters is perhaps not surprising. In the Patagonian toothfish, SOI 

798 correlates to egg and larval prevalence through its influence on ocean temperature but by the 

799 time fish recruit into the breeding population, the SOI signal has become obscured (Belchier 

800 and Collins 2008). Certainly many factors come into play between spawning and recruitment in 

801 a fish that matures 15 years later. Other, shorter-lived ground fish in the Scotia Sea have 

802 exhibited correlations to SOI, but as those authors note, these fish feed heavily on krill, a 

803 species known to track closely SOI changes to its ocean environment. Antarctic toothfish, as far 

804 as anyone knows, do not feed readily on krill; but then no one knows what toothfish eat during 

805 winter along the Ross Sea slope where krill and krill predators are abundant (Ballard et at. 

806 2011). 

807 

808 Change in toothfish prevalence 

809 Monthly Pattern in CPUE. The within-year pattern of apparent toothfish abundance (catch rate) 

810 was interesting but somewhat resists explanation given the lack of other information. The 

811 pattern shown in the 27 years of data was similar in form but with a lesser decrease in catch 
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812 rate after November than reported by Testa et al. (1985). In fact, the pattern apparent in our 

813 data was similar to one of the three years in the Testa et al. study; the other two showed a 

814 more extreme decrease for December. To explain the temporal pattern, one might be tempted 

815 to invoke some sort of migratory movement on the scale of shorebirds visiting the Arctic tundra 

816 to breed but then vacating it before food becomes unavailable, or for that matter the annual 

817 invasion and evacuation of the Ross Sea by Adelie penguins (Ballard et al. 2010). After all, 

818 McMurdo Sound tagged fish have been recaptured well to the north of McMurdo Sound, as 

819 noted (see Hanchet et al. 2008); southward movement into the southern Ross Sea by industry- 

820 tagged fish (along the slope) has been detected as well (Hanchet et al. 2010). On the other 

821 hand, toothfish could also seasonally change their position in the water column to be more 

822 vulnerable to our setline array. From personal experience with other fish species (ocean 

823 salmon), if you are not fishing at the depth where they occur you won't catch them, and their 

824 depth can change dramatically and day-to-day (D. Ainley, 100s of hours of pers. obs.). 

825 Therefore, we are inclined to explain the early spring increase in our catch rate on the basis of 

826 shifting position in the water column due to food availability. In early spring (August, 

827 September) with the sun first rising on 20 th August, there would be no reason for silverfish, and 

828 therefore toothfish, to be very high above the bottom (where our hooks are located) given a 

829 phytoplankton concentration close to zero (indicated by visibility ~80 m; Barry 1988). Any prey 

830 for silverfish would be on or near the bottom and, indeed, crystal krill are known to feed on 

831 benthic detritus during winter (Nicol et al. 2004, Deibel and Daly 2007). Therefore, perhaps all 

832 toothfish, along with their prey, are found at or near the bottom at that time, and not inclined 

833 to find our baits. A couple of months later, when the plankton blooms begin and which are 
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834 sufficiently dense to alter the 1% light level from 54 m to just a few meters deep (visibility to <6 

835 m; Barry 1988, Arrigo et al. 1998), there is ample reason for zooplankton and fish to occupy the 

836 water column. Then, the depth of toothfish tracks that of silverfish in as noted in McMurdo 

837 Sound (Fuiman et al. 2002). 

838 If there is a diminution of catch rate after November/December (perhaps in a few years; 

839 Testa et al. 1985), we hypothesize that this is related to a true decrease in toothfish abundance 

840 (especially large fish) owing to predation, which increases severely at that time as Weddell seals 

841 attempt to recover from breeding (beginning in late November), as a new cohort of seals begin 

842 to forage in the foodweb (late November), and as Ross Sea killer whales appear (early 

843 December; Ainley and Siniff 2009, Ainley et al. 2009). Dearborn (1965) and Calhaem and 

844 Christoffel (1969) reported toothfish eaten abundantly by seals as late as January in McMurdo 

845 Sound. Moreover, industrial vessels have been quite successful as well in January (Hanchet et 

846 al. 2010). Thus, the fish are present. Testa et al. (1985) showed that predation affects the 

847 spatial extent of toothfish abundance, so why not the temporal extent as well? 

848 Indeed, according to Everson (1970), the movements and annual cycles of nototheniids 

849 in the Scotia Sea appear to have evolved to avoid predation by top predators. This brings us 

850 back to the question of why it seems that as encountered by the fishery, in ice-free waters, 

851 Antarctic toothfish are found on the bottom? Sperm whales and elephant seals (summer, ice- 

852 free visitors only) are known to dive to 2400 and 3000 m, respectively, Weddell seals to 700 m, 

853 and fish-eating killer whales to 350 m (Ballard et al. 2011); of course average foraging depths 

854 are proportionately much shallower. The fishery targets the fish in waters 1000-2000 m. Living 
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855 at that depth would be out of range for at least two of these predators. In shallower waters, the 

856 smaller fish, being also cryptically colored, appear to hide within the "forest" of benthic 

857 invertebrates (Eastman and Barry 2002). Of interest, in trophic studies that looked at fish 

858 caught in deep waters along the Ross Sea slope and over sea mounts to the north, 59-64% had 

859 empty or essentially empty (food "trace") stomachs (Fenaughty et at. 2002). In contrast, surface 

860 waters over the slope are teaming with foraging predators, gorging on krill and fish (Ballard et 

861 at. 2011), and over the shelf in ice covered waters only 10% of toothfish stomachs have been 

862 found empty (Eastman 1985b). It would appear that the toothfish in the deep depths are biding 

863 their time, awaiting the return of sea ice cover and the seasonal departure of some of their 

864 predators. At this time they could be particularly vulnerable to baited hooks. 

865 

866 Decadal Pattern in CPUE. In regard to longer-term change in catch rate, no factor that 

867 we used in our analyses could explain the drop off in the CPUE of toothfish in McMurdo Sound 

868 beginning after 2001. One could hypothesize that the mega icebergs that blocked McMurdo 

869 Sound 2001-2005 were responsible (Arrigo et at. 2002), and one of the explanations favored by 

870 Hanchet et at. (2010, discussing a preliminary report of the McMurdo Sound data set), but 

871 other than a temporary effect, no other mesopredator was affected by the icebergs in a long- 

872 lasting manner (beyond 2005). In fact, Weddell seal numbers and pup production immediately 

873 returned to pre-iceberg levels (Siniff et at. 2008), Adelie penguin breeding population in the 

874 southwestern Ross Sea increased during this period and through to the present, and the 

875 prevalence of emperor penguins in McMurdo Sound increased by an order of magnitude 

876 (Landcare Research NZ, unpubl. data; Ainley pers. obs.). All of these species appear to have a 
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877 broadly similar diet during summer in waters of the southern shelf as do toothfish, principally 

878 through high consumption rates of silverfish (cf. La Mesa et al. 2004, Ballard et al. 2011). The 

879 fishing site did not change sufficiently to explain the trend either (see Hanchet et al. 2010). 

880 Coincident with the decrease in the CPUE of toothfish, toothfish length (and condition) 

881 decreased as well. We can only conclude that the industrial fishery, which targets the largest 

882 fish, is reducing their prevalence in the southwestern Ross Sea. Could it also be reducing the 

883 prevalence of "high quality" individuals (ones that grow faster and mature quickly and are in 

884 better condition), as other fisheries do when targeting the largest, oldest fish (Longhurst 2010, 

885 Ainley et al. 2011)? Apparently this is happening in the heavily fished Patagonian toothfish of 

886 the Scotia Sea, although detected over a longer time period (Shust and Kozlov 2006). Given the 

887 inability of the industrial fishery to catch many large Antarctic toothfish over the Ross Sea shelf 

888 (Hanchet et al. 2008, 2010), fishery biologists dependent on CCAMLR data would not be aware 

889 of this decreasing prevalence of large fish. Indeed, Hanchet et al. (2010) state that trends are 

890 unlikely to be found in the industrial catch data owing to the high variability in catch rates and 

891 characteristics of individual vessels. 

892 A change in the representation of large fish in the southern Ross Sea, indicative of the 

893 entire stock, would not be ecologically neutral (Longhurst 2010, Ainley et al. 2011 and 

894 references therein). How to detect this change is a challenge. Such a change was involved in the 

895 crash of Atlantic cod Gadus morhua, where the diminishment of the inshore portion of the 

896 stock occurred well before any signal was evident in the larger offshore stocks and prior to the 

897 eventual total collapse of the fishery (Longhurst 2010, and references therein). In accord with 
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898 the decreasing numbers of large fish, the Ross Sea killer whale prevalence in the southwestern 

899 Ross Sea continues to decrease (Ainley et at. 2009, unpubl. data), a pattern played out by the 

900 now endangered large-fish-eating, resident killer whales in waters of western Canada and Puget 

901 Sound (Ford et at. 2010). The killer whale social structure depends on the existence of large fish. 

902 The increase in penguin populations would be consistent as well with decreased competition 

903 with toothfish for silverfish prey. 

904 

905 Recommendations 

906 We could suggest that this toothfish fishery in the Ross Sea be closed until the natural history of 

907 Antarctic toothfish, and its silverfish prey, is much better understood, given the observation of 

908 recent declines in condition, size, and CPUE for the Antarctic toothfish of McMurdo Sound, the 

909 only monitoring conducted independent of the fishery. An easier option to implement, given 

910 that New Zealand has taken 55% of the toothfish from the Ross Sea area (Ainley et al. 2011) 

911 and that the United States imports 40% of the total world catch of "Chilean sea bass" (including 

912 Ross Sea fish; US Dept of State website), is for NZ and US fishery agencies to immediately 

913 reconstitute the scientific fishing program described herein as it was up to 2001 (annual, entire 

914 season), as well as initiate other research, e.g. monitoring of dependent species, that would be 

915 independent of the unexplained biases in the industrial sampling. Besides McMurdo Sound, 

916 vertical longlines could be deployed in an analagous fashion through the seasonally long-lasting 

917 fast ice of Terra Nova Bay as well as Edisto-Moubray bays farther north along the Ross Sea 

918 coast. Doing so would be consistent for an agency, CCAMLR (Commisison for the Conservation 
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919 of Antarctic Marine Living Resources), that has made its name as setting the pace in ecosystem- 

920 based and precautionary management and which has a formal monitoring program ( CCAMLR 

921 Ecosystem Monitoring Program; CEMP) from which management of fisheries is supposed to 

922 become effectively informed (Constable et al. 2000, Croxall and Nicol 2004, Constable 2011). 

923 Trends in the "unofficial CEMP" based on changes in numbers offish-eating killer whales, seals 

924 and penguins underway in the southern Ross Sea, by scientists engaged in pure rather than 

925 fishery-directed research, is pointing to the need for meaningful action. 
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1186 

1187 Table 1 Summary of data used to calculate changes in fish length, condition and catch per unit effort 

1188 (based on hooks deployed). If the field "total hooks fished" is blank, then no effort data were available 

1189 for that year. Other aspects of "standardizing" the effort, include: vessel = same each year (fish hut); 

1190 vessel operator = same (DeVries); fish measurer = same (DeVries); gear = same (vertical set line, number 

1191 of hooks varied, as indicated); location = same, i.e. within 2 km radius off Cape Armitage, except in a 

1192 couple of years as noted in Fig 1; depth = 415-495 m in all years. 


Year 

Dates 

24 h 
fishing 
periods 

N fish 
caught 
(measured) 

Total 

hooks 

fished 

Fish 

per 

hook 

1972 

9/15-12/13 

90 

96 (96) 



1973 

9/11-12/6 

86 

184 (184) 



1974 

9/8-12/19 

105 

470 (470) 



1975 

9/16-12/29 

105 

346 (331) 

2525 

0.137 

1976 

12/1-1/2 

33 

167 (160) 

485 

0.344 

1977 



- 



1978 

10/15-12/21 

67 

349 (341) 

1113 

0.314 

1979 

10/18-12/11 

55 

261 (255) 

938 

0.278 

1980 

10/20-12/15 

57 

202 (184) 

1373 

0.147 

1981 

10/13-12/20 

69 

268 (255) 

1986 

0.135 

1982 

10/12-12/22 

71 

499 (218) 

2939 

0.170 

1983 

10/13-12/20 

69 

553 (546) 

3155 

0.175 

1984 

10/1-12/8 

69 

205 (200) 

1338 

0.153 

1985 

10/15-12/20 

67 

318 (315) 

1690 

0.188 

1986 

11/2-12/30 

58 

151 (141) 

783 

0.193 

1987 

10/11-12/10 

60 

412 (407) 

1388 

0.297 

1988 

9/4-12/11 

98 

294 (294) 



1989 

9/7-12/6 

90 

182 (178) 

1185 

0.154 

1990 

8/31-12/8 

100 

283 (265) 

1190 

0.238 

1991 

11/19-11/30 

11 

7(7) 



1992 

9/9-11/18 

70 

32 (31) 

231 

0.139 

1993 



- 



1994 



- 



1995 



- 



1996 



- 



1997 

10/11-12/14 

64 

72(72) 

442 

0.163 

1998 

10/9-12/9 

61 

243 (243) 



1999 



- 



2000 



- 



2001 

11/14-12/11 

27 

76 (76) 

499 

0.144 

2002 

10/18-12/15 

59 

91 (91) 



2003 

10/3-12/22 

80 

17 (17) 



2004 

10/11-12/22 

74 

19 (19) 
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1193 

1194 

1195 

1196 

1197 


2005 

1 / 10 - 1/14 

5 

0 



2006 

* 

* 

* 



2007 

10 / 18 - 11/30 

73 

2 ( 0 ) 

660 

0.003 

2008 



- 



2009 

12 / 1 - 12/20 

20 

0 

117 

0.000 

2010 

10 / 26 - 12/3 

38 

10 ( 7 ) 

162 

0.043 


* Another researcher (G. Hofmann, pers. comm.) at site off Cape Armitage (same location as most others 
shown in Fig 1) in late Nov 2001, 10 24-h sets (15 hooks/set), caught 71 fish (included in the 2001 total); in 
late Nov 2006, approximately same site, 7 24-h sets (12 hooks/set), caught 0 fish. 
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1198 

1199 Table 2 Weighted regression analysis of percentile total length (A, B) and mean condition 

1200 index, K (C); shown are statistics for the two-part spline model and for the difference in slope 

1201 1972 to 1992 and from 1992 to 2010 (N = number of years analyzed). 



N 

Df 

F 

P 

Adj r 2 

A) TL 75 th percentile 






2-part spline model 

27 

2, 24 

11.07 

0.004 

0.436 

Test of difference in slope 


1, 24 

11.75 

0.002 


B) TL 95 th percentile 






2-part spline model 

27 

2, 24 

5.09 

0.014 

0.24 

Test of difference in slope 


1, 24 

7.32 

0.012 


C) Mean condition, K 






2-part spline model 

27 

2, 24 

13.77 

0.001 

0.496 

Test of difference in slope 


1, 24 

27.07 

0.0001 



1202 

1203 
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1204 Table 3 Regression analysis of total length (A, B) and condition index, K (C, D) using individuals 

1205 (N) and comparing data sets with (A, C) and without (B, D) 1991 and 1992. Shown are results of 

1206 models for 2-part spline (estimated slope up tol992 plus estimated slope sincel992, with knot 

1207 at 1992) and for the difference in the two slopes. All models included month of capture (df = 3, 

1208 for September, October, November, and December). 



N 

Df 

F 

P 

Adj R 2 

A) Total Length, all years 






full model (includes month) 

5437 

5, 5431 

30.15 

0.0001 

0.0261 

2-part spline component 


2, 5431 

41.42 

0.0001 


Test of difference in slope 


1, 5431 

25.89 

0.0001 


B) Total Length, without 1991 & 1992 






full model (includes month) 

5399 

5, 5393 

29.34 

0.0001 

0.0256 

2-part spline component 


2, 5393 

38.48 

0.0001 


Test of difference in slope 


1, 5393 

23.25 

0.0001 


C) Condition, all years 






full model (includes month) 

5403 

5, 5397 

24.66 

0.0001 

0.0214 

2-part spline component 


2, 5397 

55.94 

0.0001 


Test of difference in slope 


1, 5397 

108.61 

0.0001 


D) Condition, without 1991 & 1992 






full model (includes month) 

5365 

5, 5359 

22.98 

0.0001 

0.0201 

2-part spline component 


2, 5359 

50.66 

0.0001 


Test of difference in slope 


1, 5359 

99.46 

0.0001 
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1210 
1211 
1212 
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1213 Table 4 Summary of statistics for regression analyses with respect to fish length (75 th and 50 th 

1214 percentiles) and mean condition in relation to environmental variables, and relationship 

1215 between SOI and Max ice. N = number of years analyzed. 



N 

df 

F 

P 

Adj r 2 

TL 75 th percentile-Max ice area 4 yr lag 

17 

1, 16 

4.63 

0.048 

0.1849 

TL 50 th percentile-Max ice area 4 yr lag 

17 

1, 16 

8.78 

0.009 

0.3273 

Condition - Min ice area 8 mo lag 

20 

1, 19 

4.99 

0.039 

0.1734 

Condition - Max ice extent 3 yr lag 

18 

1, 17 

4.98 

0.040 

0.1895 

TL 75 th percentile - SOI (July-Dee) 4 yr lag 

25 

1, 24 

1.59 

0.219 

0.0242 

SOI - Max ice area same year 

29 

1, 28 

20.5 

0.0001 

0.4101 


1216 

1217 



1218 

1219 

1220 

1221 

1222 

1223 

1224 

1225 

1226 

1227 

1228 
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Figure 1 Southern McMurdo Sound showing locations where the set line was positioned 
seasonally, 1972-2010. Numeral 1 indicates position of fishing site in 2009 and 2010, as well as 
site of "Penguin Ranch", where Ponganis and Stockard (2007) recorded numerous seals 
capturing toothfish in 2003 and 2004; numeral 2 is alongside positions where fishing was 
otherwise conducted; movement of sites in part forced by shifting location of the sea ice 
runways (only the two most common runway positions shown). The dashed lines separate 
zones of fishing success attained by Testa et al. (1985; fish per day) as a function of distance 
from the major Weddell seal breeding location, which is in top right corner of map. 
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1233 



1235 Figure 3 Length-frequency of Antarctic toothfish caught in southern McMurdo Sound, 1972- 

1236 2010. Dashed line is the locally weighted regression (lowess, using running-line least squares) of 

1237 length on year. 

1238 
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1239 


1240 

1241 

1242 

1243 



percentile of length; right, change in fish length in the McMurdo Sound catch, expressed as the 
95 th percentile of length, 1972-2010. 
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1244 

1245 



1246 Figure 5 Left, mean toothfish condition and standard error among fish caught in McMurdo 

1247 Sound, monthly, 1972-2010 (n >5500 fish); right, mean catch per unit effort and standard error 

1248 by month over the period 1975-2010. 
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1250 



1251 

1252 Figure 6 Trends in mean fish condition, K, among fish caught in McMurdo Sound, 1972-2010 

1253 
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1254 



1255 Figure 7 Left, change in prevalence of axe handle fish in the McMurdo Sound toothfish catch, 

1256 1972-2010; right, relationship between proportion of axe handle fish in the catch and maximum 

1257 ice area 4 years earlier. 
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1259 



1260 

1261 Figure 8 Annual variation in fish captures (In-transformed) per 10 hooks set in McMurdo 

1262 Sound, 1972-2010. Monthly variation in capture rates, after adjusting for month of capture, as 

1263 well as best-supported model for annual variation (stair-step), are shown. 
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1265 



1266 Figure 9 Fish length as affected by maximum ice area, 4 years previous, for both 75th and 50th 

1267 percentile. 

1268 
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1269 



1271 1979-2007; right, fish condition as affected by ice extent three years previously, 1979-2007. 

1272 

1273 


72 


1274 



1275 Figure 11 Left, relationship between maximum ice area and SOI during July-December of same 

1276 year; right, trend in maximum ice extent 3 years previously during the study period (consistent 

1277 with Zwally et at. 2002, though a longer time series). 


